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Abstract 

The analytical package written in FORM presented in this paper allows the 
computation of the complete set of Feynman Rules producing the Rational terms 
of kind R2 contributing to the virtual part of NLO amplitudes in the Standard 
Model of the Electroweak interactions. Building block topologies filled by means 
of generic scalars, vectors and fermions, allowing to build these Feynman Rules in 
terms of specific elementary particles, are explicitly given in the gauge class, 
together with the automatic dressing procedure to obtain the Feynman Rules from 
them. The results in more specific gauges, like the 't Hooft Feynman one, follow as 
particular cases, in both the HV and the FDH dimensional regularization schemes. 
As a check on our formulas, the gauge independence of the total Rational contri- 
bution (Ri + R2) to renormalized S-matrix elements is verified by considering the 
specific example of the H — )• 77 decay process at 1-loop. This package can be of 
interest for people aiming at a better understanding of the nature of the Rational 
terms. It is organized in a modular way, allowing a further use of some its files 
even in different contexts. Furthermore, it can be considered as a first seed in the 
effort towards a complete automation of the process of the analytical calculation of 
the R2 effective vertices, given the Lagrangian of a generic gauge theory of particle 
interactions. 

Keywords: Electroweak interactions, NLO, Virtual Radiative Corrections, Uni- 
tarity. Rational terms 



1 Introduction 



The calculation of NLO radiative corrections to multiparticle production processes has 
achieved significant progresses in the last few years, thanks to refinements to the tra- 
ditional techniques PQ [21 [3] on the one hand, and to the introduction of new methods, 
mainly based on Unitarity [H [5] and Generalized Unitarity [6] principles, on the other. 
This has already allowed the computation of several (differential) cross-sections at NLO, 
especially for key signal and background particle scatterings and decays of interest at 
colliders, whose signatures can potentially through light on the mechanism underlying 
the Electroweak (EW) Symmetry Breaking process [UIB]. 

A generic 1-loop amplitude can be decomposed as a linear combination of known scalar 
integrals, with up to 4 external legs, plus a residual Rational part R. In the framework of 
the approaches based on the Unitarity of the S-matrix, the first piece is Cut-Constructible 
(CC), i.e. it can be reconstructed by properly cutting the loop amplitude in tree- level 
like sub-amplitudes, thus conceptually reducing the complexity of the calculation of 1 - 
loop integrals to a simpler computation of tree-level diagrams. On the other hand, when 
working in approaches based on 4 integer dimensions, the R terms, or at least some of 
them as clarified in the following, cannot indeed be calculated just in terms of tree- level 
diagrams, but a full 1-loop computation is needed. This can be performed according 
to the traditional Feynman diagram approach [3 [ID]- An alternative strategy consists 
in making use of Unitarity-based on-shell recursion relations [H] , leading to a recursive 
evaluation of R by means of bootstrapping techniques [12]. On the other hand, it is 
worth observing that the R terms can be put on the same footing as the CC part at 
the price of introducing a larger number of integer dimensions, as worked out in the d- 
dimensional extensions of the Unitarity approach [13] . Despite the attractive elegance of 
the d- dimensional Unitarity formulations, in this paper we choose to follow a variation 
of the first approach, thus avoiding to introduce extra integer dimensions (and, as a 
consequence, to extend to these higher dimensions the external particle wave functions), 
and just allowing ci = 4 + e dimensions in the dimensional regularization scheme we use 
to regularize the divergencies appearing during the computation of 1-loop integrals (even 
in case of finite tensor integrals, singularities may arise in the tensor reduction process 
and need to be regularized). We worked in the framework of the OPP method [H], one 
of the Unitarity inspired algebraic and universal (i.e. independent from the model of 
particle interactions at hand) procedures first introduced to automatically calculate the 
CC part of any 1-loop amplitude. In the OPP approach, the R terms can be organized in 
two classes, according to their origin: the Ri terms, arising from the mismatch between 
the 4-dimensional part of the numerator and the c?- dimensional denominator including 
the poles of the propagators in the regularized integrand of 1-loop integrals, and the R2 
terms, arising instead from the e-dimensional part of the numerator. Both these classes 
are thus a residual effect of the dimensional regularization scheme one introduces on the 
integrands to calculate 1-loop integrals. In the framework of OPP, it has been shown that 
the Ri terms are closely related to the CC part, and can indeed be obtained numerically 
at the same time of the last one [1^ . Unfortunately, however, this property does not apply 
to the R2 terms, that, instead, need a dedicated computation. Anyway, the problem of 
the calculation of the R2 part of a generic 1-loop amplitude can be solved by observing 
that only 1-particle irreducible diagrams with up to 4 external legs can contribute to R2, 



1 



due to the exclusively ultraviolet nature of these terms, as proven in [TO]. In other words, 
R2 terms will never contribute to the 1/e and pole parts of 1-loop amplitudes, whose 
nature is instead completely infrared. Thus, taking into account that the number of the 
contributing diagrams is limited, it is possible to calculate, once and for all for the theory 
of interaction at hand, all possible R2 effective vertices, up to 4-external legs, and then 
use these effective vertices as building blocks for computing the R2 contribution to each 
specific 1-loop amplitude [15]. Given the set of all external particles (with their momenta 
and their quantum numbers) identifying a 1-loop helicity amplitude, it will be enough 
to consider all tree-level diagrams joining them, which include one (and only one) R2 
effective vertex. This recipe has been adopted in the HELAC-l-loop code [TB], where 
the R2 effective vertices entering QCD 1-loop corrections [T7] have been implemented, 
and recursion relations have been adopted to evaluate the corresponding contributions to 
the amplitudes. Furthermore, our R2 Feynman Rules, together with (off-shell) recursive 
relations for tensor integrals, have been applied to the numerical computation of the R2 
contribution to 1-loop multigluon amplitudes in a tensor reduction framework [T5] . 

Throughought this paper we denote by topology a set of lines (propagators) connect- 
ing a set of points (vertices), by generic diagram a topology filled by means of generic 
Standard Model (SM) scalar s, vector v and/or fermion / fields, and by particle diagram 
a generic diagram whose components have been further specified in terms of selected SM 
particles (like H, Z, e~ , etc.). After a formal definition of the R terms in Section H] 
and our comments on the gauge choice in Section [31 we present our analytical results 
for the R2 contributions to all leading generic 1-loop 1-particle irreducible diagrams with 
up to 4 external legs arising in the SM of EW interactions, diagram by diagram, and 
the FORM code we have written to obtain from them the R2 effective vertices involving 
real specific particles, to be used in actual calculation of 1-loop amplitudes, in Section HI 
These analytical results are of interest for people aiming at better understanding the 
nature and the properties of the R2 terms. We include a test of the reliability of our 
code together with considerations concerning the gauge invariance of R in Section [5l 
and we draw our conclusions in Section |6l Additional information about the notation 
used in our code is available in the Appendix and in the README file within the package. 



2 The R contribution to 1-loop amplitudes 

In this Section we formally define the R terms, and in particular the R2 class, briefly 
reviewing how they appear in the computation of 1-loop amplitudes. Our starting point 
is the integrand of a generic m-point amplitude, written as 

A{q) = „ /^^l, , with A = (g + - , (1) 

where q is the loop momentum, q + Pi and mj are the momentum and the mass of the 
i-th loop particle {i = 0, l,....m — 1). 

A dimensional regularization procedure is then introduced, at the aim of regularizing 
the singularities in the integrand amplitude, in order to be allowed to evaluate its integral. 
In renormalizable gauge theories, the number of ultraviolet (UV) divergent integrals is 



2 



finite, and tlie UV divergencies are re-absorbed, after integral calculation, in renormalized 
quantities. The infrared (IR) singularities, instead, cancel when combining together the 
real and the virtual contribution to the amplitude at each fixed order in the perturbative 
expansion in terms of the coupling constant, according to the KLN theorem Thus 
the results for 1-loop amplitudes will, in general, include a finite part plus a divergent 
part, of IR origin only. 

We choose to work in c? = 4 + e dimensions (dim). The extension to rf-dim of the 
integrand of the amplitude, N{q) — )■ N{q), is achieved through the transformations: 

Qf^^q,, = q,i + qii, (2) 
7m ^ 7m = 7m + 7a ' (3) 
Qfiu 9 flu = Qf^u + QfLv , (4) 

where we have denoted quantities in d-dim with a bar, and we have explicitly shown the 
e-dim part of each quantity, including tildes. The quantities in e-dim are orthogonal with 
respect to the quantities in 4-dim, thus, in developing the expression of N{q), one has to 
take into account relations like 

qfiv^ = qX^ (5) 

'gf.ug^" = 0, (6) 

and similar ones. 

One can thus rewrite the numerator in d-dim as a sum of a 4-dim part plus a residual 
part. The R2 terms correspond to the integral of the e-dim part of N{q), i.e. 

1677^ J D0D1D2 Dm-l 

To understand the origin of the CC and the Ri terms, one can instead expand the 
4-dim part of the numerator, N{q), in terms of 4-dim denominators Di, according to the 
universal OPP decomposition [H], and then observe that d-dim denominators instead 
of 4-dim denominators appear in the dimensionally regularized integrand of the ampli- 
tude [20]. By rewriting these d-dim denominators in terms of the 4-dim denominators of 
the decomposition according to 

1 A q'\ 1 
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one can see that the first part of the expression in the parentheses lead to the CC part 
of the amplitude (all terms in 4-dim, as tree- level diagrams), whereas the second part of 
the expression in the parantheses causes the appearence of the Ri contribution to the 
amplitude: 

i?i = ^^-.^OT^ / . (9) 
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3 Choice of a gauge 



One of the crucial points to actually perform the computation of the R2 terms is the 
choice of a convenient gauge. We work in the class of gauges. More precisely, we 
consider the generalized gauges, characterized by a three-parameter ^a, ^z, ^ gauge 
fixing term. The standard R^ gauges can be obtained as the particular case corresponding 



The main difference between the gauges in this class and the unitary one, is the ap- 
pearence in the first ones of unphysical scalars and Fadeev-Popov-DeWitt ghosts, as loop 
particles. The unitary (even called physical) gauge instead, comes out as the particular 
limit C,, 00 and C,a finite number (these two independent limits are possible 

because the first one does not fix the electromagnetic gauge invariance), corresponding 
to unphysical particles becoming extremely heavy and decoupling from the theory^ In 
fact, the unitary gauge condition eliminates 3 of the 4 degrees of freedom in the scalar 
Higgs doublet, re-emerging as longitudinal spin states of the Z and gauge bosons 
acquiring a mass. Thus, one has to deal with a reduced number of particles, i.e. with the 
physical ones only. This reduced number of particles is the reason why, when working 
at tree-level, unitary gauge is often preferred. However, although loop calculations can 
indeed be performed even in this gauge, this is not convenient, due to the fact that the 
unitary gauge is not manifestly renormalizable. In particular, the UV behavior of the 
theory is not properly manifest and appears worse than it really is. The reason is that 
the expression of the propagators of the massive gauge bosons in the unitary gauge 



does not fall to zero for large momenta g — )■ 00. 

On the other hand, i?^ gauges are covariant gauges, where the unphysical degrees of 
freedom in the Higgs doublet do not disappear. However, the gauge fixing part of the 
Lagrangian eliminates the possibility of their mixing with gauge bosons, thus allowing the 
development of a perturbation theory, and assigns them a mass which is related to the 
one of the gauge bosons by the conditions = y/^niz, = ^JJjriw- The propagators 
of the vector bosons in the generalized R^ gauges have the expression 



that, in the limit ^, 00, tends to eq. ( [TU]) valid in the unitary gauge. Eq. (ITT]) 

however, when compared to eq. (ITU]) , has the advantage of being applicable both to 
massive and to massless vector bosons. Furthermore, in the limit q — t- 00, it goes to zero, 
thus the UV behaviour of the theory is not spoiled by the choice of this gauge. Actually, 
we worked both in the generalized i?^ gauges and in the unitary gauge, and we explicitly 
verified that the R2 effective vertices obtained in the latter have expressions far more 
complicated [22j, as expected on the basis of the previous considerations!! 

^ Actually, the Fadeev-Popov ghosts do not completely decouple from the theory. As ^, — )■ 00, 
there are some surviving pieces (Lee- Yang terms) coming from the H-ghost-ghost vertices proportional 



^When working in the Unitary gauge, at the purpose of calculating the R2 (and the Ri) effective 



to u = iz = i- 




(10) 




(11) 
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4 Procedure followed in the calculation of the R2 ef- 
fective vertices and structure of the code 

We worked in the SM of the EW interactions, and we calculated in an analytical way, by 
using FORM [23], the general R2 contributions corresponding to all possible 1-particle 
irreducible graphs, with up to 4 external legs, written in terms of generic scalars s, vectors 
V and fermions / (as explained in the Introduction, we denote the graphs where generic 
fields, like s, v and /, appear instead of specific particles, like if, Z, e~, etc., as generic 
diagrams throughout this paper). 

The starting point for the evaluation of the R2 contribution to each generic diagram 
is the corresponding R2 integrand expression (under the integral sign in eq. ([7])), built as 
suggested in Section [2|, once fixed a gauge, on the basis of the considerations in Section [31 
Feynman rules written in terms of generic s, v and / according to Denner conventions [21] , 
enter the numerator of this integrand. The integration can then be performed analitically 
by using conventional techniques, involving Feynman parametrization followed by Wick 
rotation. The basic formulas for the angular and radial integrations are given by |H] 



The integration always leads to the disappearence of all divergencies involving negative 
powers of e. One can thus safely take the limit e — )■ just after it. 

By following this procedure, we obtain the R2 analytical expressions corresponding 
to the generic diagrams listed in the xxxxgentop.h files present in the package, with 
xxxx = ss, vs, vv, ff , sff , vff , sss, vss, svv, vvv, ssss, ssvv and vvvv. The name 
of each of these files refers to the nature of the external particles and the results are 
presented generic diagram by generic diagram (e.g. the file vssgentop.h includes all 
R2 contributions to the 3-point functions including one vector and two scalars, generic 
diagram by generic diagram, where each generic diagram differs from the others according 
to the topology and/or to the nature of the fields internal to the loop). 

Only the generic diagrams giving rise to an R2 contribution different from zero are 
listed. The main criterion for establishing if a diagram can contribute or not to R2 is 
power counting in the loop momentum. To every integral of the type 



we associated an index d' = l + s + l — m. Integrals with ci' > can give a contribution, 
while integrals with rf' < vanish [25]. By checking in each diagram the number of 
available loop momenta in the numerator and by counting the number of denominators, 
we could say beforehand if they may contribute or not. All gauges in the class, included 
the 't Hooft Feynman one, have the same power counting properties. On the other hand, 

vertices, we took the limit ^, — ^ cx) and — !■ 1 before the integration, according to the observations 
in 




(12) 



(13) 




(14) 
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in the unitary gauge the situation is different and more diagrams can contribute than in 
the previous case, at least as for power counting, since the piece of the propagator ffTOj) 
proportional to q^q^ just includes two more powers of momentum in the numerator and 
thus passes more easily the power counting test. This is counterbalanced by the fact 
that, due to the reduced number of degrees of freedom, many diagrams are absent (e.g. 
all those including a ssv vertex, since two physical s don't couple with a v and diagrams 
including unphysical Goldstone s are absent). Another criterion is the existence or not 
of terms proportional to e or (f' in N{q). In the case of ghosts for example, the loop 
momenta in the ghost-ghost-vector Feynman rules have an index that comes from an 
external vector and therefore is 4 dimensional. When such an object is contracted with 
a loop momentum, it cannot produce q^ terms and thus can not contribute to R2 termsH 

The generic diagrams giving a non-null R2 contribution in the i?^ gauges are shown 
in Fig. [1], [2] and [3|, where each of them is presented in a single topological configuration, 
just for compactness. In fact, only selected topologies have been explicitly dressed to 
generate generic diagrams, i.e. we explicitly considered only one topology associated to 
each fixed number of internal and external legs, and then different topologies (and the 
corresponding generic diagrams) have been obtained from the first ones by proper non- 
cyclic permutations of the momenta of the external particles, at a fixed configuration of 
the internal particles. In case of external bosons, besides the momenta, even the Lorentz 
indices have to be properly permuted, to obtain one topology/ diagram from another. The 
schemes showing our convention corresponding to the topologies explicitly considered for 
2- leg (bubbles and tadpoles), 3- leg (triangles and bubbles) and 4- leg (boxes, triangles and 
bubbles) diagrams, are shown in Fig. |H 

Each R2 contribution included in a xxxxgentop . h file is labelled by EFFVERyyyyLzzzz, 
where yyyy correspond to the ordered list of the generic external particles (i.e. the first 

particle corresponds to the pi momentum, the second one to p2, etc ), L stays for "loop" , 

and zzzz is the ordered list of the generic particles running inside the loop, where the 
first internal particle is, by convention, the one joining external particle 1 with external 
particle 2, and so on, up to the last internal particle joining the last external particle to 
the first external one. 

Furthermore, the R2 contribution corresponding to each generic diagram is presented 
as an explicit FORM function of the external and internal particles. As shown in Fig. HJ 
we associated two labels to every internal particle, each of them in correspondence with 
the different vertex to which it is joined. The reason of this choice is better accomo- 
dating charged particles and fermions (an internal fermionic line corresponds to both 
an antifermion emerging from a vertex and to a fermion entering in another vertex), by 
taking into account that in each vertex all momenta are supposed to be incoming. This 
notation seems a little bit redundant, however we found it very useful in making clearer 
the code. 

Furthermore, in case of fermionic loops, just for practical purposes, we explicitly 
distinguish two generic diagrams corresponding to the same topology, obtained one from 
the other at a fixed configuration of external legs, by simply changing the direction of 
the fermion flow (clockwise in the ff , iff, ffff diagrams, and anticlockwise in the 
ffac, fffac, ffffac diagrams). In case of the ssvv diagrams, an additional couple of 

Ghost loops do not contribute to R2 terms. These loops however may contribute to the Ri part of 
the amplitude, where have to be taken properly into account. 
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Figure 1: Non null contributions to the ss, vs, vv and ff R2 effective vertices in the 
generalized gauges, with generic finite ^, ^z, C,a- The corresponding analytical for- 
mulas associated to selected topologies of each generic diagram are included in the files 
xxxxgentop .h, with xxxx = ss, vs, vv and ff, respectively. 

topologies (clockwise and anticlockwise) appears, corresponding to the exchange of the 
second scalar with the first vector, that has to be considered separately (see Fig. |3]). 

The symmetry factors have not been included at this level of the calculation. They 
have instead been restored at the following step, i.e. when replacing generic fields with 
specific particles, accomplished by means of do-loop procedures as explained below. 

All global variables, functions and constants used in the package are declared in 
the file variables. h, including some explanations of their meaning. In particular all 
SM particles, with their mass, charge and isospin, are listed (see also Section |H] on our 
notations). 

The expressions of all EW vertices, taken from Ref. [2^, are encoded in a specific 
file, named f illvertices .h, organized in procedures. Each procedure dresses (i.e. fills) 
a particular generic vertex (e.g. svv) with real specific particles. Special care has been 
taken in the treatment of the vss, vvv and vvvv vertices. In fact, these vertices are 
not symmetric under the exchange of two generic particles, so one has to distinguish 
different possibilities. In case of generic diagrams involving vvvv vertices, we explicitly 
obtain different R2 formulas according to the nature (charged W or neutral V) of the four 
vectors entering the vertex, and to their position in the diagram as external or internal 
particles. In general, one has to distinguish 4 cases corresponding to the same topology, 
as shown in Fig. O the case of an internal {W^, W~) and an external [W^, W~) couple, 
the one of an internal {W, V) and an external {W, V) couple, the one of an internal 
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Figure 2: Non null contributions to the sff , vff , sss, vss, svv and vvv R2 effective 
vertices in the generalized gauges, with generic finite ^, ^z, ^a- The corresponding 
analytical formulas associated to selected topologies of each generic diagram are included 
in the files xxxxgentop.h, with xxxx = sff, vff, sss, vss, svv and vvv, respectively. 
For all diagrams including a vvvv vertex see the comments in Fig. O 

{V, V) and an external {W~^, W~) couple, and the one of an internal {W~^, W~) and an 
external {V, V) couple. In practice, the first two cases and the second two cases can be 
grouped together, corresponding to the same formulas. We explicitly distinguish these 
possibilities in the formulas presented in the files. A further configuration has also to be 
taken into account in building the W^W~W^W~ R2 effective vertex, corresponding to 
an internal {W~^, W~^) and an external {W~ , W~) couple or viceversa. 

The procedure to fill generic masses with the ones of real particles and some procedures 
that relate particles with their charge conjugated, just for internal use, are also included in 
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Figure 3: Non null contributions to the ssss, ssvv and vvvv R2 effective vertices in 
the generalized gauges, with generic finite ^, ^z, ^a- The corresponding analytical 
formulas associated to selected topologies of each generic diagram are included in the 
files xxxxgentop.h, with xxxx = ssss, ssvv and vvvv, respectively. For all diagrams 
including at least a vvvv vertex see the comments in Fig. O 

the f illvertices .h file. All procedures included in the f illvertices .h file are valid in 
the 't Hooft Feynman gauge. The additional multiplicative ^, factors appearing in some 
of the expressions of particle masses when considering the generalized gauges instead 
of the 't Hooft Feynman one, are always already included directly into the R2 formulas 
associated to the generic diagrams listed in the xxxxgentop.h files. In particular, we 
introduced a csif function, which gives the ^ parameter associated to each particle under 
consideration. To control its values the user is allowed to modify the procedure f illcsi 
in the f illvertices .h file. Further warnings about the validity and the applicability of 
some of these formulas are explicitly mentioned in the files. 

The generic diagrams are then dressed by means of specific particles, in such a way to 
sum together all 1-loop contributions corresponding to the same real external particles, 
by considering all different non-equivalent topologies. The R2 effective vertices are built 
this way, by using the do-loop procedures included in the xxxxf eynrules . f rm files, with 
xxxx= ss, vs, vv, ff , sff , vff , sss, vss, svv, vvv, ssss, ssvv, vvvv. 

These do-loop procedures have the advantage of leading to all results for each class of 
effective vertices joining specific particles in an automatic way, without the need of con- 
sidering each effective vertex as a particular case, to be treated separately from the others 
according to the specific external particles attached, and of explicitly drawing all Feyn- 
man particle diagrams contributing to it. The most external do- loops run on the external 
particles, whereas the other ones run on the internal particles. All possible combinations 
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Figure 4: 2-, 3- and 4-point 1-particle irreducible topologies explicitly considered for the 
calculation of the R2 contributions included in the xxxxgentop.h files. The symbols 
Cj (extjfla in our files), with j = 1,2,3,4, denote external particles. Two indices are 
associated to each internal particle ik (k = 1,2,3,4), corresponding to the two vertices to 
which it is connected. So ejik (ejintfcf la in our files) denotes the k-th internal particle 
incoming in the vertex in which also the j-th external particle is entering. According to the 
notation used to write down the Feynman rules, all momenta in all vertices are supposed 
to be incoming. Different topologies have been obtained from these ones by non-cyclic 
permutations of the external particles, together with their momenta and their Lorentz 
indices (omitted for simplicity in this figure), if present, for each fixed configuration of 
the internal ones. 



of internal and external particles are considered at the do-loop level, then, the ones that 
do not exist in nature, are simply discarded due to the lack of the corresponding vertices 
in the f illvertices .h procedures. The do-loop procedures have been tested against a 
completely independent filling procedure, that works effective vertex by effective vertex, 
merely consisting in hand writing down, one by one, all diagrams contributing to each 
effective vertex, on the basis of the output of a FeynArts [26] computation, and summing 
the corresponding R2 contributions together. 

By simply changing the values of some control variables, the user has the possibility 
to modify the output, e.g. to study the impact of different contributions to each effective 
vertex joining specific particles, by selecting the box, triangle or bubble diagram contribu- 
tions only, or the fermionic loop contributions only, or to choose among different gauges 
in the class, in particular the 't Hooft Feynman and the Landau gauge, or among 
different dimensional regularization schemes, i.e. the Four Dimensional Helicity (FDH) 
scheme, or the 't Hooft Veltman (HV) one (different schemes differ for the treatment of 
the polarization vectors of the particles: in the FDH scheme all helicities are treated in 
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Figure 5: Diagrams in terms of generic fields, including a vvvv vertex, considered sepa- 
rately in our work. The charged vector bosons are denoted by whereas the neutral 
vector bosons are denoted by V. The first two diagrams are grouped together since they 
give rise to the same R2 formulas in the xxvvgentop.h files, as well as the second two. 

4 dimensions, in the HV scheme unobserved polarization vectors are continued to d dim, 
whereas the observed ones are kept in 4 dim [27]). 

The particular cases here mentioned can be obtained by simply changing the values of 
specific variables included at the end of each do-loop procedure, as suggested in Table [H 
which shows the already implemented options. 

The fermionic contribution is gauge independent, i.e. it is the same in all gauges, and 
this is the reason why we allow its selection separately from the other contributions. The 
part of the results proportional to Xhv is also gauge independent, i.e. the choice of a 
regularization scheme can be safely performed in an independent way with respect to the 
gauge choice. 

The package is available as a zipped tar archive at the web address 
http://www.ugr.es/~garzelli/R2SM. After decompressing it, all files with a .frm ex- 
tension can be runned by using a working installation of FORM [23], downloadable from 
the web at the address http://www.nikhef.nl/~forin, that needs to be preinstalled by 
the user. In particular, we extensively used the version 3.2 (May 2008) of the FORM 
package. We verified that the same output is produced when running the code with 
FORM version 3.3 (August 2010), now available on the web. 
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choice of the dimensional regularization scheme 


parameter 


option 


lambdahv generic 
lambdahv = 
lambdahv = 1 


default 
FDH scheme 
HV scheme 



choice of the gauge 


parameter 


option 


csia, csiz, csi generic 
csia = csiz = csi 
csia = csiz = csi = 1 


generahzed gauges, default 

1-parameter gauges 
't Hooft-Feynman gauge 



selection of different R2 contributions 


parameter 


option 


dummyb 


dummy c 


dummyd 


dummyf 


1 


1 


1 


1 


default, all contributions (boxes, triangles, bubbles) 
to each effective vertex are kept 


1 











only bubble contributions are selected 
(excluding bubble fermionic loops, if any) 





1 








only triangle contributions are selected 
(excluding fermionic triangles, if any) 








1 





only boxes contributions are selected 
(excluding fermionic boxes, if any) 











1 


only fermionic-loop contributions are selected 


1 


1 


1 





only non-fermionic-loop contributions are selected 



Table 1: Options already implemented for the control variables the user is allowed to 
tune at the end of the do-loop procedures in the xxxxf eynrul . f rm files, to select a 
particular gauge and regularization scheme, or to enlighten particular contributions to 
the i?2 effective vertices coming from specific set of diagrams. 



5 A check of the gauge invariance of the R contri- 
bution to (renormalized) S-matrix elements: the 

H — )^ 77 decay 

The Green functions in general depend on the gauge choice, but renormalized S-matrix 
elements must be gauge independent, since they correspond to physical observable quan- 
tities. 

The contributions to renormalized S-matrix elements of the CC terms and of the R 
terms are separately gauge invariant. As a particular check of this point and of the 

correctness of our calculations, we have explicitly proven that the total R contribution, 
Ri + R2, to the S-matrix element corresponding to the iJ — > 77 physical decay process at 
1-loop is gauge independent. The code we have written to produce this result is included 
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in the subdirectory RtotHAA of our package. All new variables specific to Ri calculations, 
are included in the additionalvariables .h file. The non trivial Ri contributions to 
1-loop generic diagrams including as external particles a generic scalar and two generic 
vectors, are presented in the svvRlgentop.h file. The calculation of the Ri terms has 
been performed in a similar way to the one of the R2 terms, the main difference being the 
fact that only the 4-dim part of the numerator N{q) can contribute to Ri. The Ri part 
of the tensor integrals involved has been extracted by applying the Passarino-Veltman 
reduction technique [1] in a straightforward way, by disregarding the contribution of the 
scalar functions A, B, C, D which, by definition, already contribute to the CC part of 
the amplitude. Ri effective vertices corresponding to specific svv configurations involving 
physical external particles, can be obtained by dressing the generic diagrams mentioned 
above with specific particles, thanks to a do-loop procedure, and by conveniently summing 
the corresponding contributions together, in a way analogous to the one followed to 
build the R2 effective vertices. The do-loop procedure we have used to calculate the 
R = Ri + R2 effective vertices for all specific svv combinations is presented in the 
file svvRf eynrules . f rm. Finally, in the file haa.frm, we show that the R contribution 
to the S-matrix element of the H'y'y decay process is gauge invariant, i.e. the result is 
independent of any gauge parameter (and it remains the same even in the limit ^, ^ 00 
and 1: we have repeated this calculation in the unitary gauge, making these limits 

before the integration, and obtained the same final result for the R part of the S-matrix 
element, also in agreement with Ref. |28]). 

6 Conclusions 

We have presented an analytical package, written in FORM and available on the web, 
to calculate the R2 effective Feynman rules by means of which one can compute the R2 
contribution to 1-loop amplitudes for whichever process in the SM of the EW interactions. 
As recommended in the literature, we have chosen to work in the class of gauges, 
identified by finite values of the ^, ^z, parameters, due to the good suitability of this 
class to loop calculations, since the UV behaviour of the theory is not spoiled by the 
particular expression of the propagators of the massive gauge bosons in this class. We 
have explicitly verified that this general recommendation also apply to the calculation of 
the R2 contribution alone, to 1-loop amplitudes: even if, from the technical point of view, 
analytical calculations of the R part of the amplitude can be performed as well in the 
unitary gauge, the gauges have by far better properties in terms of the compactness 
and the simplicity of the final R2 analytical formulas with respect to the unitary gauge. 

We have also offered a check of the gauge-invariance of the total R contribution, 
Ri + R2, to the renormalized S-matrix element for the if 77 process. 

The package is modular, and different pieces can be conveniently reintegrated in other 
calculations. The user is allowed to play with some parameters, in order to specify a 
particular gauge in the R^ class or a particular dimensional regularization scheme, and to 
enlighten different partial contribution to the R2 effective vertices. These options allow 
everybody to reproduce, as particular cases, some of the results presented in our previous 
papers [221 EH] , in a straightforward way. 

Finally, we think that our effort can be considered a first seed towards the more 
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complex task of automatically recovering, given a Lagrangian for whichever model of 
particle interactions, all R2 effective Feynman Rules for the theory at hand. 
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8 Appendix: Notation 

We closely follow the conventions in the paper [23]. We adopt the Bjorken-Drell or 
"mostly minus" metric convention 77^,^ = diag{+l, —1, —1, —1). 

We consider 1-loop amplitudes with at most 4 external legs. The momenta of the 
external particles are denoted by pi, p2, Ps, Pi and are all supposed to be incoming. In 
the 2 point effective vertices, the Lorentz indices associated to the external bosons whose 
momenta are pi and P2{= —pi), are respectively denoted by a and 6. In the 3 point 
effective vertices, the Lorentz indices associated to the external bosons whose momenta 
are pi(= —p2 — p3), p2 and ps, are denoted by r, u and respectively. Finally, in the 
4 point effective vertices, the Lorentz indices associated to the external bosons whose 
momenta are pi, p2, Ps and P4,{= —pi — p2 — p3) are denoted by a, (3, t and x- 

We worked in the generalized gauge, by allowing 3 different parameters ^a, and 
^ (csia, csiz and csi in our files) in the gauge fixing term (usually taken to be equal 
in the standard gauge). The \hv (lambdahv in our files) dependence, identifying the 
dimensional regularization scheme, has been made explicit in all results. We checked that 
the contribution proportional to Xhv for each effective vertex is always invariant with 
respect to gauge transformations. 

In our FORM files, we labelled the particles of the EW sector of the SM according 
to the notation of Table [2l Actually, we limited to consider only one weak-isospin family 
of quarks and one weak-isospin family of leptons. Anyway, the results can be extended 
to three families in a straightforward way, by just working on the fermionic part of the 
do-loop procedures in all files with a .frm extension. The number of quark colors can 
be fixed by the user by a proper choice of the value of the variable ncol at the end of 
each do-loop procedure. The contribution of fermionic loops is independent of the gauge 
choice, since all propagators and coupling constants involving fermions are the same in 
all gauges. Thus, for the parts of our final results including fermion loops, and for their 
generalization to three weak-isospin families, the reader can safely refer to the analytical 
formulas already presented in our previous paper [29] . 

The mass of each SM particle is denoted by adding an "m" before the name of the 
particle (e.g. mhsca is the mass of the Higgs scalar, mwvec is the mass of the W gauge 
boson). The pole masses of the unphysical scalars x and (j)^ are reduced to the masses 
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Particle content of the EW sector of the SM 


generic particle 


specific particles 


notation in our FORM files 


scalars s 




hsca, chisca, phiplussca, phimensca 


vectors v 


A, Z, W+, W 


avec, zvec, wplusvec, wmenvec 


ghosts fp 


Ua, Uz, M^y, 

Ua, Uz, u^, Uw 


gavec, gzvec, gwplusvec, gwmenvec 
gavecbar, gzvecbar, gwplusvecbar , gwmenvecbar 


charged leptons / 


e, /i, r 
e, /i, f 


ele, muele, tauele 
elebar, muelebar, tauelebar 


neutral leptons / 




nuele , numuele , nutauele 
nuelebar, numuelebar, nutauelebar 


up-type quarks / 


u, c, t 
u c, i 


uquark, cquark, tquark 
ubarquark , cbarquark , tbarquark 


down- type quarks / 


d, s, b 
d s,b 


dquark , squark , bquark 
dbarquark, sbarquark, bbarquark 



Table 2: Notation used in this work for the particle content of the EW sector of the 
SM. The Faddeev-Popov-De Witt ghosts are only relevant for the calculation of the Ri 
terms, whereas all other particles enter in the calculation of both the R2 and the Ri 
contributions to 1-loop amplitudes. 

of the corresponding massive gauge bosons by = y/^mz and = yjjrriw We 
also widely apply the relation mz = mw/cos9w, and thus express our results in terms 
of mw- The sine and cosine of the Weinberg angle are denoted by sinwei and coswei, 
respectively, and the e coupling constant is denoted by eem. 

The third component of the weak-isospin of each fermion is denoted by including the 
prefix "13" before the name of the particle: iSnuele (= 1/2) corresponds to the Uf,, 
iSele (= -1/2) corresponds to the e~, i3u (= 1/2) and i3d (= -1/2) correspond to the 
u and d quark quantum numbers, respectively. 

The projector operators fi+ = and n_ = are denoted by omegaplus and 
omegaminus, respectively. As for the treatment of 75, the reader can refer to the comments 
already presented in Ref. PU] . 
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